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Abstract 

Background: Back pain represents a substantial burden globally, ranking first in a recent assessment among causes 
of years lived with disability. Though back pain is widely studied among working age adults, there are gaps with 
respect to basic descriptive epidemiology among seniors, especially in the United States. Our goal was to describe 
how pain, function and health-related quality of life vary by demographic and geographic factors among seniors 
presenting to primary care providers with new episodes of care for back pain. 

Methods: We examined baseline data from the Back pain Outcomes using Longitudinal Data (BOLD) registry, the 
largest inception cohort to date of seniors presenting to a primary care provider for back pain. The sample included 
5,239 patients > 65 years old with a new primary care visit for back pain at three integrated health systems 
(Northern California Kaiser-Permanente, Henry Ford Health System [Detroit], and Harvard Vanguard Medical 
Associates [Boston]). We examined differences in patient characteristics across healthcare sites and associations of 
patient sociodemographic and clinical characteristics with baseline patient-reported measures of pain, function, and 
health-related quality of life. 

Results: Patients differed across sites in demographic and other characteristics. The Detroit site had more African- 
American patients (50%) compared with the other sites (7-8%). The Boston site had more college graduates (68%) 
compared with Detroit (20%). Female sex, lower educational status, African-American race, and older age were 
associated with worse functional disability as measured by the Roland-Morris Disability Questionnaire. Except for 
age, these factors were also associated with worse pain. 

Conclusions: Baseline pain and functional impairment varied substantially with a number of factors in the BOLD 
cohort. Healthcare site was an important factor. After controlling for healthcare site, lower education, female sex, 
African-American race, and older age were associated with worse physical disability and all of these factors except 
age were associated with worse pain. 

Trial registration: Clinical Trials.gov NCT01 776242; Registration date: June 13, 2012. 



* Correspondence: jarvikj@uw.edu 

'Department of Radiology, University of Washington, 325 Ninth Ave., Box 
359728, Seattle, WA 98104-2499, USA 

department of Neurological Surgery, University of Washington, Seattle, WA, 
USA 

Full list of author information is available at the end of the article 

O© 2014 Jarvik et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
BlOlVlGCl C^ntrBl Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain 
Dedication waiver (httpy/creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, 
unless otherwise stated. 



Jarvik et al. BMC Musculoskeletal Disorders 201 4, 1 5:1 34 
http://www.biomedcentral.com/1471 -2474/1 5/1 34 



Page 2 of 1 1 



Background 

Back pain is the most common reason worldwide for years 
living with disability (YLDs) [1]. In the United States in 
2010, low back pain was the top contributor to YLDs, out- 
ranking diseases such as chronic obstructive pulmonary 
disease and diabetes [2]. Although back pain is common 
among adults aged 65 years and older ("seniors"), it has 
been under-studied in this age group [3,4]. Recently, sev- 
eral publications from groups outside the U.S. have begun 
to explore this issue [5-7]. Nonetheless, relatively little is 
known concerning how back pain, function, and health- 
related quality of life (HRQoL) differ according to patient 
sociodemographic and clinical characteristics, especially 
among the elderly in the United States. 

In this study, we addressed these knowledge gaps 
using baseline data from the Back pain Outcomes using 
Longitudinal Data (BOLD) cohort of patients aged 
65 years and older initiating a new episode of care for 
back pain [8]. Our primary focus was to determine the 
associations of patient age, sex, race, ethnicity, and edu- 
cation with patient-reported measures of pain, function, 
and HRQoL. A unique strength of BOLD is the geo- 
graphic diversity of the cohort participants. Therefore, 
we evaluated the magnitude of the differences in the 
characteristics of back pain patients across recruitment 
sites, the strength of the association between patient 
characteristics and baseline measures, and whether asso- 
ciations between patient characteristics and patient- 
reported outcomes are consistent across the three re- 
cruitment sites. 

Methods 

Institutional review board (IRB) approval 

The study was approved by the Institutional Review Boards 
(IRBs) of all the participating institutions (University of 
Washington, Harvard Vanguard, Henry Ford Health 
Systems, Northern California Kaiser-Permanente). 

Patients and setting 

We previously described the details of the registry [8]. 
In brief, using healthcare system electronic databases, 
we prospectively identified patients aged 65 years or 
older with a primary care visit for back pain in the prior 
3 weeks at 3 integrated healthcare systems: Harvard 
Vanguard (Boston), Henry Ford Health Systems (Detroit), 
and Kaiser-Permanente Northern California. We excluded 
patients who had visits for back pain in the 6 months prior 
to the index visit. Our other primary inclusion and ex- 
clusion criteria were the International Classification of 
Diseases, Ninth Revision, Clinical Modification (ICD-9- 
CM) diagnostic codes [9] listed in Additional file 1. We 
approached patients in-person, by telephone or mail 
and performed and administered all baseline measures 
either in-person or by telephone within 3 weeks of their 



index visit, verifying eligibility and obtaining consent 
at that time. We enrolled patients from March 2011 
through March 2013. 

Measures 

We obtained the following patient-reported measures at 
baseline: 

Sociodemographic characteristics: age, sex, race, ethni- 
city, education, employment status, marital status, smok- 
ing status and whether they had a lawyer involved with a 
back-related claim. 

Pain-related characteristics: a) duration of current epi- 
sode of back/leg pain (less than 1 month, 1-3 months, 
3-6 months, 6-12 months, 1-5 years, more than 5 years); 
b) average back pain intensity and average leg pain inten- 
sity in the past week on 0-10 numerical rating scales 
where 0 = 'no pain' and 10 = 'pain as bad as you can im- 
agine'; c) the 24-item Roland-Morris Disability Question- 
naire (RMDQ) [10], modified to specify disability related 
to either back or leg pain; and d) Brief Pain Inventory 
(BPI) Activity Interference Scale [11,12] (consisting of 7 
0-10 ratings of how much pain interferes with general ac- 
tivity, mood, ability to walk, normal work, relations with 
other people, sleep and enjoyment of life). 

Psychological distress: the 4-item PHQ-4 (0-12) meas- 
ure of anxiety and depressive symptoms has been demon- 
strated to be a general marker of psychological distress 
[13]. PHQ-4 scores of 6 or greater have been recom- 
mended as "yellow flags" and scores of 9 or greater as "red 
flags" for presence of a depressive or anxiety disorder [14]. 

HRQoL: EuroQol-5D (EQ-5D), a preference-weighted, 
quality of life index (0-1) consisting of five dimensions 
(mobility, self-care, usual activities, pain/discomfort, and 
anxiety/depression) and a visual analog scale of current 
HRQoL [15]. 

Falls: Number of falls in the past 3 weeks and how 
many resulted in injury, from the Behavioral Risk Factor 
Surveillance System (BRFSS) survey [16]. 

Recovery Expectations: Patients rated their confidence 
that their back and/or leg pain would be completely gone 
or much better in 3 months on a scale from 0 = 'not at all 
confident' to 10 = 'extremely confident' [17,18]. 

Statistical analysis 

Due to the large sample size of the registry and known 
geographic and sociodemographic differences across the 
three recruitment sites, statistical comparisons of base- 
line variables across sites (ANOVA for continuous vari- 
ables and Chi-square for categorical variables) yielded 
statistically significant differences at p = 0.05 on all mea- 
sures except for age. We use boxplots (mean +/- stand- 
ard deviation) and descriptive statistics to characterize 
RMDQ and back pain scores by patient sociodemo- 
graphic characteristics and ICD-9-CM diagnostic codes 
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assigned at the index visit. To examine the independent as- 
sociations of baseline demographics with patient-reported 
outcome measures, we created separate multivariable lin- 
ear regression models for RMDQ and back pain scores, 
adjusting for all available baseline demographic variables 
as well as recruitment site. We report unstandardized 
model coefficients and two-sided p-values without adjust- 
ment for multiple testing, because our research goal was 
to characterize the healthcare site differences rather than 
to test specific a priori hypotheses. 

Results 

Patient enrollment 

Figure 1 depicts the flow of patients in the study. Of 
13,376 patients identified as potentially eligible, we were 
unable to contact 15%, 15% were ineligible, and 27% de- 
clined to participate or to complete the baseline ques- 
tionnaire. The remaining 5,239 patients (39% of patients 
identified for screening) enrolled and completed the 
baseline study measures. The mean (SD) number of days 
between the index visit and the baseline assessment was 
14.6 (5.3). 



Patient differences across sites 

Although patients from the three sites were similar in 
age, they differed on all other sociodemographic vari- 
ables (Table 1). Notably, the Detroit site was character- 
ized by the highest proportions of African-Americans 
(50%), females (70%), and current smokers (13%). The 
Boston site had more highly educated patients (68% with 
college degree), and more patients who were married 
(80%), followed in both cases by the Northern California 
site (39% with college degree, 58% married). 

Detroit patients also differed from those at the other 
sites with respect to many of the patient-reported out- 
come measures (Table 2): they were less likely to have 
had back pain for less than 3 months, they rated their 
back and leg pain as more intense, they reported 
greater pain interference with activities, and they re- 
ported lower HRQoL, as measured by the EQ-5D. On 
the PHQ-4 measure of psychological distress [13,14], 
the mean score at the Boston site was lower than 
those at the other sites, but mean scores at each site 
were well below the "yellow flag" cutpoint of 6 points. 
The mean score on the recovery expectation rating 
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13,376 patients identified for screening through ICD-9 code 
2,037 Ineligible (>3 weeks from index visit) 
505 Not reached 
3,499 Declined further screening 
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1,973 Patients ineligible for one or more reasons: 
666 Incomplete screening 
403 Prior visit for back pain in last 6 months 
399 Had a major medical comorbidity 
372 Did not speak English 
25 1 Had prior back surgery 
211 Recent visit not back-related 

62 Failed cognitive screening assessment 

61 Other 
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5,382 Patients were eligible 

133 Patients were not included 

74 did not consent for further participation 
59 did not complete the baseline questionnaire 



Figure 1 Consort diagram depicting patient flow from initial screening through enrollment. 
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Table 1 BOLD cohort sociodemographic characteristics overall and by healthcare system 


Demographic variable 


Combined 
N = 5239* 


N. CA 
N = 3164 


Detroit 
N = 967 


Boston 
N = 1108 


Age, mean (SD) years 


73.8 ±6.9 


73.7 ±6.8 


73.8 ± 6.9 


73.8 ±7.1 


Sex, N (%) 










Male 


1 852 (35%) 


1 1 54 (36%) 


289 (30%) 


409 (37%) 


Female 


3387 (65%) 


2010 (64%) 


678 (70%) 


699 (63%) 


Hispanic, N (%) 


305 (6%) 


270 (9%) 


21 (2%) 


14 (1%) 


Race, N (%) 










Black or African American 


797 (1 5%) 


226 (7%) 


485 (50%) 


86 (8%) 


Native American Indian/Alaskan/Hawaiian or other Pacific Islander 


44 (1%) 


40 (1%) 


4 (0%) 


0 (0%) 


Asian 


1 98 (4%) 


171 (5%) 


1 8 (2%) 


9 (1 %) 


Caucasian 


3843 (73%) 


2430 (77%) 


432 (45%) 


981 (89%) 


Other 


299 (6%) 


246 (8%) 


26 (3%) 


27 (2%) 


Education, N (%) 










Less than high school graduate 


314 (6%) 


144 (5%) 


154 (16%) 


16 (1%) 


High school graduate or obtained a GED/Vocational, technical/trade school 


1450 (28%) 


794 (25%) 


382 (40%) 


274 (25%) 


Some college 


1 275 (24%) 


978 (31%) 


234 (24%) 


63 (6%) 


Four year college graduate 


1 290 (25%) 


605 (19%) 


89 (9%) 


596 (54%) 


Professional or graduate degree 


894 (1 7%) 


633 (20%) 


105 (11%) 


156 (14%) 


Marital, N (%) 










Married/living with a partner 


3191 (61%) 


1 850 (58%) 


451 (47%) 


890 (80%) 


Separated/divorced 


602 (11%) 


423 (13%) 


147 (15%) 


32 (3%) 


Never married and presently single 


256 (5%) 


1 63 (5%) 


69 (7%) 


24 (2%) 


Widowed 


1 1 75 (22%) 


716 (23%) 


299 (31%) 


160 (14%) 


Employment, N (%) 










Working full-time/part-time 


602 (11%) 


367 (12%) 


166 (15%) 


69 (7%) 


Retired (not due to ill health) 


4250 (81%) 


2535 (80%) 


913 (82%) 


802 (83%) 


Retired or disabled because of ill health 


146 (3%) 


71 (2%) 


12 (1%) 


63 (7%) 


Other 


220 (4%) 


1 77 (6%) 


11 (1%) 


32 (3%) 


Legal representation for back pain, N (%) 










Yes 


31 (1%) 


23 (1%) 


8 (1%) 


0 (0%) 


No 


5196 (99%) 


3136 (99%) 


956 (99%) 


1104 (100%) 


Smoking status, N (%) 










Never Smoked 


2901 (55%) 


1570 (50%) 


675 (70%) 


656 (59%) 


Quit smoking over a year ago 


1 999 (38%) 


1422 (45%) 


169 (17%) 


408 (37%) 


Current smoker/quit less than a year ago 


324 (6%) 


1 64 (5%) 


121 (13%) 


39 (4%) 


Days between index visit and baseline assessment, mean (SD) 


14.6 ±5.3 


16.0 ±4.6 


16.3 ±5.3 


8.9 ±2.7 



•Numbers vary slightly for items due to slight variation in response rates with overall missing items being less than 19 



was somewhat lower at the Detroit site than at the 
other sites. 

Relationship of patient sociodemographic characteristics 
to patient-reported pain, physical disability, and HRQoL 

Figure 2 shows RMDQ and back pain intensity scores at 
each site stratified by patient age, race, sex, and educa- 
tion. RMDQ scores increased with age, from a mean 



(SD) of 9.2 (6.6) among those aged 65-69 to 10.7 (6.1) 
for those older than 85 across sites. The oldest age 
group (>85) had worse RMDQ scores at all three sites. 
The average pain duration also increased with age, with 
31% of patients in the youngest age group (65-69) hav- 
ing had pain more than one year versus 44% in the old- 
est age group (>85). The EQ-5D index was similar for all 
age groups [mean (SD) = 0.76 (0.17)]. 
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Table 2 Patient-reported measures overall and by healthcare system site 



Baseline symptoms and outcomes 


Combined 


N. CA 


Detroit 


Boston 




N = 5239* 


N = 3164 


N = 967 


N = 1108 


RMDQ (0-24), mean (SD) 


9.5 ± 6.4 


9.7 ± 6.0 


12.8 ±5.9 


6.4 ± 6.4 


Back pain duration, N (%) 










< 1 month 


1 749 (33%) 


1215 (38%) 


230 (24%) 


304 (27%) 


1 - 3 months 


1015 (19%) 


715 (23%) 


137 (14%) 


163 (15%) 


3 - 6 months 


344 (7%) 


243 (8%) 


46 (5%) 


55 (5%) 


6-12 months 


313 (6%) 


1 65 (5%) 


96 (10%) 


52 (5%) 


1 - 5 years 


776 (1 5%) 


399 (1 3%) 


1 97 (20%) 


180 (16%) 


> 5 years 


1037 (20%) 


426 (1 3%) 


261 (27%) 


350 (32%) 


Back pain intensity (0-10), mean (SD) 


5.0 ±2.8 


4.9 ±2.7 


6.2 ± 2.8 


4.5 ± 2.7 


Leg pain intensity (0-10), mean (SD) 


3.4 ±3.3 


3.5 ±3.2 


4.2 ±3.6 


2.6 ±3.1 


BPI Interference (0-10), mean (SD) 


3.3 ±2.5 


3.4 ± 2.5 


3.8 ±2.6 


2.5 ± 2.3 


PHQ-4 (0-12), mean (SD) 


1.6 ±2.5 


1.9 ±2.5 


1.8 ±2.8 


0.4 ± 1 .6 


EQ-5D index (-0.1-1), mean (SD) 


0.76 ±0.1 7 


0.76 ±0.1 7 


0.69 ±0.20 


0.81 ±0.16 


EQ5D VAS (0-100), mean (SD) 


74.3 ± 1 8.4 


73.6 ± 18.9 


74.0 ±18.4 


76.6 ± 1 6.9 


Falls in past 3 weeks, N (%) 










Patients with one or more fall 


385 (7%) 


288 (9%) 


81 (8%) 


16 (1%) 


How many of these falls caused one or more injury**? 


1 79 (46%) 


121 (42%) 


51 (63%) 


7 (44%) 


Back Pain Recovery Expectations (0-10), mean (SD) 


5.5 ±3.7 


5.7 ±3.6 


4.9 ±3.7 


5.1 ±4.0 


• *Sample sizes vary slightly for items due to slight variation in response 


rates with overall 


missing items being less than 0.3°/ 


b. Complete data were available for 



RMDQ, back pain intensity, and BPI interference. 

• **The number and percent of patients in whom a fall caused at least one injury where injury was defined as limiting regular activities for at least a day or 
requiring a visit to a doctor. 



African-Americans reported worse pain and pain- 
related physical disability than did Caucasians at all 
study sites (Figure 2). The mean (SD) RMDQ score was 
12.1 (6.5) among African-Americans compared with 
8.9 (6.3) for Caucasians. The proportion of African- 
American patients whose pain duration was longer than 
a year was 42% compared with 34% for Caucasians, and 
the mean scores for back and leg pain were 0.9 and 1.0 
points worse, respectively, for African-Americans than 
for Caucasians after adjusting for site. 

We observed these relationships at each site as well as 
for the overall cohort. For example, although African- 
Americans at the Detroit site reported greater physical 
disability than did African-Americans at the Northern 
California site (RMDQ mean (SD) = 13.6 (5.9) versusl0.6 
(6.2)) or Boston (7.3 (7.1)), within each site, RMDQ scores 
were higher for African-Americans than for Caucasians 
(Detroit 13.6 (5.9) vs. 12.0 (5.9), Northern California 10.6 
(6.2) vs. 9.4 (6.0), and Boston 7.3 (7.1) vs. 6.2 (6.3)). 

Patients from different racial groups differed in educa- 
tion. More Caucasians (47%) than African-Americans 
(19%) were college graduates. 

Because of the small number of Hispanic patients at 
Detroit (21) and Boston (14), we report data from His- 
panics from Northern California only. The mean RMDQ 
score was 1.2 points higher (worse) for Hispanics than 



for non-Hispanic Caucasians [10.7 (6.0) vs. 9.5 (6.0)]. 
Hispanics also reported greater back pain intensity, ac- 
tivity interference, and depression/ anxiety. 

Women reported significantly worse function, pain, 
and activity interference (Figure 2e and f). However, 
men and women did not differ in back pain duration, 
PHQ-4 scores, EQ5D scores, or recovery expectations. 

Although there were marked differences across sites in 
levels of education (Figure 2g and h), both across and 
within sites, patients with less than high school educa- 
tion reported the greatest pain duration, back pain, leg 
pain, and physical disability. Among patients with less 
than high school education, 51% reported pain duration 
greater than 1 year, compared with only 27% among 
those with a college degree. 

Patients who had retired or were disabled due to ill 
health had worse pain and physical function than pa- 
tients still working. Patients who had retired not due to 
ill health also had somewhat worse pain and physical 
function than did patients still working. 

Association of diagnosis with pain, function, and HRQoL 

Table 3 lists the six most commonly recorded diagnosis 
codes (see Additional file 1 for a complete list of ICD-9- 
CM Diagnosis Codes). Lumbago was used more than 
three times as frequently as the next most commonly 
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Figure 2 Roland-Morris Disability Questionnaire (RMDQ) and average back pain intensity in past week by age, race, sex and education 
stratified by site, a, c, e, and g display the baseline measure of back-related physical disability (the Roland-Morris disability questionnaire) and 
b, d, f, h display the baseline measure of pain (a numerical rating scale of average back pain in past week), stratified by site and compared with 
key demographic variables, a and b depict baseline measures by age; c and d depict baseline measures by race; e and f depict baseline measures 
by sex; and g and h depict baseline measures by education. Each outcome by demographic boxplot displays the within-group mean (horizontal 
in the center of the box) +/- the within-group standard deviation (upper and lower edges of the box). The vertical lines display the within-group 
range of scores. 



Table 3 ICD-9-CM Diagnosis codes 



Diagnosis codes 


N (%) 


RMDQ 




Back pain intensity 


Leg pain intensity 


Non-Specific: lumbago, backache, sprains and strains 


3560 (67.95%) 


4.0 9.0 1 


5.0 (9.5 ± 6.5)* 


3.0 5.0 7.0 (5.1 ±2.8) 


0.0 2.0 6.0 (3.0 ± 3.3) 


Back and Leg: back pain with radiation, sciatica 


1091 (20.82%) 


5.0 10.0 


15.0 (10.0 ±6.2) 


3.0 5.0 7.0 (4.9 ± 2.9) 


3.0 5.0 7.0 (5.1 ± 3.0) 


Spinal Stenosis 


288 (5.50%) 


5.0 10.0 


15.2 (10.0 ±6.4) 


3.0 5.0 6.0 (4.6 ± 2.6) 


0.0 5.0 6.0 (4.0 ± 3.2) 


Others 


300 (5.73%) 


2.0 7.0 1 


3.0 (8.0 ± 6.6) 


3.0 5.0 7.0 (4.7 ± 2.7) 


0.0 0.0 4.0 (2.2 ± 2.9) 



*25 t percentile, median, 75* percentile (mean ± standard deviation). 
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used code. Patients with the diagnosis codes that indi- 
cated leg involvement (sciatica, back pain with radiation, 
and stenosis) had higher leg pain intensity ratings than 
did patients with diagnoses that did not include leg in- 
volvement (lumbago, backache, strain and sprain). Pa- 
tients with diagnosis codes indicating back pain with leg 
involvement as well as stenosis diagnoses codes had the 
highest RMDQ scores (mean (SD) = 10.0 (6.4) and 10.0 
(6.2), respectively) (Figure 3). 

Multivariable analysis of association of baseline 
characteristics with RMDQ and pain 

The multivariable analysis demonstrated strong as- 
sociations of site, employment status, education and 
duration of symptoms with RMDQ scores after ad- 
justment for other baseline demographic variables 
(Table 4). We also observed modest associations of 
race, age, gender and smoking status with RMDQ 
scores. We observed similar patterns of association for 
back pain intensity ratings, although, unlike for RMDQ 
scores, back pain intensity was negatively associated 
with the number of days since the index visit and was 
not associated with age. Neither marital status nor His- 
panic ethnicity was associated with RMDQ scores or 
back pain intensity. 



Discussion 

This is the largest cohort to date of adults aged 65 and 
older with a new episode of care for back pain. Our re- 
sults show that there are important differences in pain 
intensity, physical disability, and health-related quality of 
life across different healthcare sites and across different 
patient age, sex, and racial subgroups. These findings are 
of importance to both researchers and clinicians. For re- 
searchers comparing the effectiveness of interventions in 
observational studies, our findings emphasize the im- 
portance of adjusting for any patient baseline differences 
between treatment groups in factors that are associated 
with outcome measures of pain, function, and HRQoL. 
For example, education level appears to be important in 
this regard, but often is not reported or adjusted for in 
such studies. Our findings also indicate that identifiable 
subgroups of older patients may differ substantially at 
the time of initiating a new episode of care for back pain. 
This suggests the potential value of applying different inter- 
ventions tailored to these different subgroups. Moreover, 
the baseline characteristics of patients and the healthcare 
system environment in which they are located should be 
considered when evaluating treatment outcomes. 

We observed sizable differences in baseline patient- 
reported measures across recruitment sites, with patients 
from the Detroit site worse on most measures as 




Primary ICD9 
^| NorvSpecific 
\—\ Back and Leg 
r~ H Spinal Stenosis 



\ CA 



Detroit 



Boston 



Figure 3 RMDQ and back pain intensity by ICD-9-CM code stratified by site, a displays the baseline measure of back-related physical 
disability (the Roland-Morris disability questionnaire) and b displays the baseline measure of pain (a numerical rating scale of average back pain in 
past week), stratified by site and compared with primary diagnosis code for the index visit using ICD-9-CM. Each outcome by demographic 
boxplot displays the within-group mean (horizontal in the center of the box) +/- the within-group standard deviation (upper and lower edges of 
the box). The vertical lines display the within-group range of scores. 
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Table 4 Multivariable linear regression for dependent variables of baseline RMDQ and back pain NRS on patient 
characteristics 



Independent 

i/ariaHl^ 

VO I Id DIC 


Levels 


Dependent variable: RMDQ (0-24) 


Dependent variable: 
back pain intensity 
(0-10) 








Coef* (95% CI) 


P-value 


Coef* 


(95% CI) 


P-value 


Site 


Boston 


Ref 


<0.001 


Ref 




<0.001 




Detroit 


D.22 (4.00, D.oo) 




1 .26 (0.96, 1 .56) 






N. CA 


oJU U.//, o.o/) 




0.59 (0.35, 0.83) 




Age 


(Per year) 


r\ r\ i / r\ r\r\ r\ r\r\ 

0.03 (0.00, 0.05) 


0.048 


0.00 (- 


■0.01, 0.01) 


0.983 


Sex 


Male 


Ref 


<0.001 


Ref 




<0.001 




Female 


1.10 (0.73, 1.46) 




0.50 (0.34, 0.67) 




Hispanic 


Yes 


Ref 


0.578 


Ref 




0.087 




No 


-0.24 (-1.10, 0.62) 




-0.34 


;-0.73, 0.05) 




Race 


Black or African American 


Ref 


0.011 


Ref 




<0.001 




Native American Indian/Alaskan/Hawaiian or other 
Pacific Islander 


-0.13 (-1.93, 1.67) 




-0.05 


;-0.88, 0.77) 






Asian 


— U./y {— I 0.6oj 




-0.50 


;-0.94, -0.06) 






Caucasian 


— O./o (—1.2/, — 0.2b) 




-0.72 


;-0.96, -0.49) 






Other 


o.io (— o./y, i . i o) 




-0.51 


;-0.94, -0.07) 




Education 


Less than high school graduate 


Oaf 


<U.UU 1 


Ref 




<0.001 




High school graduate/obtained a GED/Vocational, 
technical, or trade school 


110 / 10-7 mm 

— 1 .1 3 (—1 .8/, —0.39) 




-0.35 


;-0.69, -0.02) 






Some college 


-2.20 (-2.96, -1.44) 




-0.40 


;-0.75, -0.06) 






Four year college graduate 


-2.99 (-3.78, -2.20) 




-0.80 


[-1.16, -0.44) 






Professional or graduate degree 


-3.04 (-3.85, -2.23) 




-1.05 


;-1.42, -0.68) 




Marital 


Married/living with a partner 


Ref 


0.386 


Ref 




0.086 




Separated/divorced 


0.20 (-0.34, 0.73) 




0.07 (- 


■0.17, 0.32) 






Never married and presently single 


0.45 (-0.31, 1.21) 




0.25 (- 


■0.10, 0.60) 






Widowed 


0.33 (-0.12, 0.78) 




0.25 (0.04, 0.45) 




Employment 


Working full-time/part-time 


Ref 


<0.001 


Ref 




0.003 




Retired (not due to ill health) 


0.86 (0.33, 1 .38) 




0.30 (0.06, 0.54) 






Retired or disabled because of ill health 


3.90 (2.82, 4.99) 




0.87 (C 


.38, 1.37) 






Other 


0.60 (-0.33, 1 .54) 




0.08 (- 


■0.34, 0.51) 




Legal rep. 


Yes 


Ref 


0.059 


Ref 




0.115 




No 


-2.00 (-4.07, 0.08) 




-0.76 


(-1.71,0.19) 




Smoking status 


Never Smoked 


Ref 


<0.001 


Ref 




0.012 




Quit smoking over a year ago 


0.39 (0.04, 0.74) 




0.13 (- 


■0.03, 0.29) 






Current smoker/quit less than a year ago 


1.46 (0.77, 2.15) 




0.45 (0.14, 0.77) 




ICD9 code 


Non-specific 


Ref 


<0.001 


Ref 




0.051 




Back and Leg 


1 .00 (0.59, 1 .42) 




-0.04 


;-0.23, 0.15) 






Spinal Stenosis 


1.08 (0.35, 1.81) 




-0.46 


;-0.79, -0.13) 






Others 


-0.75 (-1.46, -0.05) 




-0.14 


;-0.47, 0.18) 
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Table 4 Multivariable linear regression for dependent variables of baseline RMDQ and back pain NRS on patient 
characteristics (Continued) 



Pain duration 0-1 month 
1 - 3 month 
3-6 month 
6-12 month 
1 - 5 years 
> 5 years 

Time from index visit to assessment (per day, 0-21) 



Ref 

0.62 (0.16, 1.08) 
0.61 (-0.07, 1.30) 
1 .07 (0.35, 1 .78) 
1.80 (1.29, 2.31) 
2.03 (1.56, 2.50) 
-0.04 (-0.08, -0.00) 



<0.001 



0.029 



Ref 

0.13 (-0.08, 0.34) 
0.33 (0.02, 0.65) 
0.42 (0.10, 0.75) 
0.46 (0.23, 0.69) 
0.53 (0.31, 0.74) 
-0.07 (-0.08, -0.05) 



<0.001 



<0.001 



*Multivariable regression model coefficients are on the scale of the dependent variable and are not standardized. 



compared to the Boston and Northern California sites. 
Site differences in disability and pain persisted even after 
controlling for demographic factors available to us. It is 
possible that unmeasured socioeconomic differences are 
responsible at least in part for the site differences. Detroit 
was suffering from a severe economic recession during 
the study time-frame, and depressed economic conditions 
have been shown to be associated with poorer well-being 
and quality of life in the elderly [19]. Our findings 
emphasize the limitations of single-site observational 
studies, which may not be generalizable to other settings. 

Our finding that less educated patients reported worse 
function is concordant with a review by Dionne and col- 
leagues of education and back pain [20]. They speculated 
that education may be a marker for other factors, such as 
ability to adapt to stress, access to healthcare, occupational 
factors, and behavioral/environmental factors. All patients 
in our study had access to healthcare. However, it is pos- 
sible that less educated patients also were financially disad- 
vantaged and may have delayed care due to concerns 
about having to pay co-insurance costs. Less educated pa- 
tients may be economically disadvantaged and under more 
psychosocial stress, which could affect health outcomes. 
Another possibility is that the association of lower educa- 
tion with worse health-related measures reflects the cumu- 
lative effects of social disadvantage on disease burden [21]. 

The finding that women present with worse pain and 
physical disability than men is consistent with prior re- 
search in largely younger populations [22-24]. 

While there was increasing pain-related physical disabil- 
ity (RMDQ score) with age, pain severity was not clearly 
associated with age. Prior studies in conditions other than 
back pain have yielded conflicting findings regarding the 
association between age and pain severity. Creamer found 
no association between age and pain severity in a sample 
of patients with knee osteoarthritis [25]. Thomas [26] and 
Parsons [27] found increasing pain severity with increas- 
ing age for a variety of musculoskeletal conditions. In a 
nationally representative sample of Medicare beneficiaries, 
pain reporting did not vary by age [28]. 

Somewhat surprising was that, independent of age, pa- 
tients who had retired for reasons other than ill health 



also had worse pain and physical disability than those 
patients still working. While this might be due in part to 
the healthy worker effect [29], other authors have found 
retirement associated with a variety of symptoms such 
as declining mobility and daily activities, and declining 
mental health [30], including depression [31]. 

In our multivariable analysis, we found that even after 
adjusting for site, education, pain duration and other 
factors, African-American race was associated with 
worse baseline physical disability and pain at presenta- 
tion for back pain-related care. This finding is consistent 
with other studies that found worse pain-related disabil- 
ity in cohorts of African-Americans compared with 
other races [32-35]. Thus, our observation of worse pain 
and physical disability among African-Americans com- 
pared with Caucasians could be explained by different 
coping strategies, or could be a result of residual con- 
founding in our multivariable analysis. 

Patients with diagnosis codes indicating leg involve- 
ment or spinal stenosis reported slightly worse physical 
disability. This is concordant with other studies indicat- 
ing that patients with leg involvement have more severe 
pain and physical disability than those without [36-38]. 

Adjusting for other variables, as compared with patients 
who never smoked, those who were former smokers re- 
ported somewhat higher levels of physical disability and 
pain, and current smokers reported even higher levels of 
disability and pain. Substantial prior research has linked 
smoking to worse back pain outcomes, and one study of 
patients with spine-related back or leg pain found that 
compared with patients who had never smoked, current 
smokers reported greater pain; in longitudinal analyses, 
compared with patients who continued to smoke, those 
who quit reported significantly greater improvement in 
pain [39]. 

It is worth emphasizing that our large sample size 
allowed us to detect what were frequently relatively 
small differences in patient reported measures between 
subgroups. The importance of the magnitude of these 
differences on an individual level is uncertain. 

One of the limitations of our project is that we 
enrolled nearly two-thirds of the cohort from a single 
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site - Northern California Kaiser-Permanente. A second 
limitation is that nearly all of the Hispanics were from one 
site (Northern California) and a majority of African- 
Americans were from another site (Detroit). This distribu- 
tion of patients limits the overall generalizability of our 
findings. Differences observed in self-reported measures 
and outcomes may reflect site-specific differences that are 
based on local healthcare system or patient-specific fac- 
tors. Because patients at the 3 sites had different sociode- 
mographic characteristics, we will need to control for 
these factors in future analyses. A third limitation is that 
we enrolled 39% of patients initially identified as poten- 
tially eligible. Because these patients did not complete 
questionnaires, we cannot further characterize this non- 
enrolled group, but we acknowledge that they potentially 
limit generalizability. 

Another limitation is that this study was not designed to 
determine sociodemographic differences between sub- 
groups with low back pain, but rather was designed to 
examine the natural history of back pain among seniors at 
three integrated health systems. Subgroup differences may 
reflect institutional and other local factors that we did not 
measure, such as income and co-pays, that could influence 
access to healthcare and hence utilization and outcomes. 

While ours is the first cohort of seniors with back pain 
assembled from a primary care setting in the United 
States, there have been similar cohorts assembled inter- 
nationally. The Back Complaints in the Elders (BACE) 
[40] group is a consortium of investigators from the 
Netherlands, Australia and Brazil who are assembling 
similar but smaller cohorts, planned to be around 750 
patients per national cohort. The first of these to be 
published was the Dutch BACE cohort that enrolled 675 
patients [7]. Their inclusion criteria were similar to ours, 
recruiting primary care patients with a new episode of 
care for back pain. They included slightly younger pa- 
tients (>55 years old). They also had a slightly shorter 
allowed interval between the index visit and when they 
contacted patients, allowing a maximum of 2 weeks 
compared with 3 weeks for our study. These studies will 
provide an opportunity to compare and contrast the 
presentation, diagnosis and treatments of seniors with 
back pain between the U.S. and other countries. 

Conclusions 

In summary, we have presented the BOLD baseline re- 
sults describing the demographic and baseline patient- 
reported measures. Site was strongly associated with 
baseline patient reported measures. It is unclear whether 
site differences were due to geographical factors, socio- 
cultural factors specific to a particular area, or factors 
specific to the healthcare system that oversees the care 
of these patients. Adjusting for site reveals continued as- 
sociation among key demographic variables of race, age, 



and education with presenting pain. Our results lay the 
groundwork for future studies involving this cohort. 

Additional file 



Additional file 1: Baseline diagnostic categories and ICD-9 CM 
diagnosis codes included in each. The diagnoses in italics accounted 
for approximately 80% of subjects. 
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